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ABSTRACT 


Tupaz et al. (1978) formulated a model to study the be- 
havior of waves in an unstable jet which varied downstream. 
The linearized barotropic vorticity equation was solved 
numerically. Waves were forced with a fixed period on the 
eastern boundary and a radiation condition was applied on 
the western boundary. 

In this thesis, cyclic boundary conditions are used on 
the eastern and western boundaries. The numerical solutions 
show amplitude growth in time which is approximately expo- 
nential. The solutions are normalized and the wave struc- 
ture 1S obtained during three time segments. The wave 
Structures are not the same during these segments, but they 
are very similar in areas of large amplitudes. In these 
regions the behavior is very similar to that obtained by 


Tupaz et al. (1978). 
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I. INTRODUCTION 


A jet, regardless of its latitudinal position in the at- 
mosphere, may contain regions of large vorticity gradients 
where the necessary condition for barotropic instability is 
locally satisfied. Synoptic-scale moving disturbances occur- 
ring at the level of a moderately strong easterly jet south 
of the Tibetan high near 200 mb level appear to arise from 
barotropic instability of the mean flow (Krishnamurti, 197la, 
1971b). If the observed disturbances are the result of baro- 
tropic instability, they will extract energy from the mean flow. 

Many investigators have studied the linear stability of 
barotropic zonal flows over the years. Krishnamurti's obser- 
vational and Colton's (1973) numerical studies both showed 
that the zonal variation of the jet apparently had significant 
effects on the dynamic behavior of the transient disturbances. 

Tupaz et al. (1978) further extended the study of the 
zonal variation of the basic flow. A numerical model was 
developed based on the linearized non-divergent barotropic 
vorticity equation on a beta-plane. This study examined the 
dynamics of transient barotropic waves in a region of variable 
mean wind. The mean zonal wind was an easterly hyperbolic 
Secant-squared, or Bickley, jet and the mean meridional wind 
was derived so that the mean flow would be non-divergent. 
Using the numerical model, a periodic forcing was applied on 


the inflow (eastern) boundary to generate waves which were 
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allowed to propagate through the channel and out through the 
western boundary. As the waves moved, they reacted to the 
local stability properties of the mean flow, whereas at each 
point the fields varied periodically. The results, which 
were obtained from long term integration, were compared with 
the parallel flow theory by constructing a simple mechanistic 
analytical model which incorporated the local stability con- 
cept of the parallel flow theory. 

In this thesis, the model developed by Tupaz et al. (1978) 
is modified by changing the boundary conditions. The quanti- 
ties on the eastern and western ends of the channel are set 
equal to enforce cyclic continuity. This allows a wider 
variety of wind profiles to be considered as long as they 
are periodic in x. The use of these boundary conditions will 
also improve our understanding of the mathematical behavior 


of waves inamean flow with downstream variation. 
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II. BASIC EQUATIONS 


A. GENERAL FORMULATION 


This study employs the barotropic vorticity equation to 


describe motion in a periodic channel. 


The development which 


follows Tupaz et al. (1978) begins with the barotropic vor- 


ticity equation on a beta (8) plane: 


(22) 


= + ou = + Vv $2 + 6V > 0. =D 
where 
Here, 8 is the constant value of the north-south gradient of 


Coriolis parameter computed at 10° latitude. 


O- “iS>a~ forerngd 


function which represents the non-barotropic and diabatic 


effects which are required to maintain the basic state vor- 


ticity field, and De 


1s a frictional coefficient. 


The mo- 


tion is assumed to be non-divergent so that a stream function 


(W) can be introduced which is defined by 


The vorticity becomes 


UN 
tt 

<J 
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The non-divergent barotropic vorticity equation may now be 


written, 
2 2 
ee ay avew , av acy re 2 
Be a= ay ee oy ee (2.5) 


The linearization of eq. (2.5) is accomplished by separat- 
ing the stream function into its mean (bar) and perturbation 


(prime) quantities of the form: 


Wika yo) e Uke) Peo Key te (2.6) 


where |w'l < < |w{ . Substituting (2.6) into (2.5) and 


dropping the quadratic terms in w' leads to 


pee? Je 


Q7] a 
<j! 
Q 
J 
= 
ale 


_ ce ee ee 
st YY OX aX Oy Vay dy aX au 
ay' 8 25 4g BW 2p gy 
+S ay Vb + Bay DV (2.7) 
It is assumed that Q' = 0 and Q is defined in such a way 


that wp is a steady-state solution of (2.5). Equation (2.7) 


can also be rewritten in Jacobian form such that 


2 ) a im 2 l i j 2 re) 
V + a isi al (¥.V p ) rz el (y »V W) = 6 st 
or 
~ DEV W 3 (2.8) 
where 1 (a,b) = oa ob ‘os da dd . 


Equation (2.8) is a Poisson equation for the tendency of the 


J 
perturbation stream function (Se). 
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Equation (2.7) is solved in a channel defined by 


i 
a 
[A 
« 


< D 


Xy SX S Xe 


The boundary conditions of the north-south walls y = + D 


are 
yp' = 0. (2.9) 


Cyclic boundary conditions are used in the east-west direc- 
tion. Tupaz et al. (1978) applied a time-periodic forcing 

on the eastern boundary to introduce waves into the region 

from the east, and used a radiation condition for both ten- 
dency and vorticity at the western boundary to allow the 


waves to move smoothly through the boundary. 


B. ZONALLY VARYING BASIC FLOW 
The basic velocity field is defined by an easterly Bick- 


ley jet 


u(x,y) = U(x) sech* (a4) -U,=- 


Here, d(x) is a characteristic length scale of the jet, U(x) 


(2.10) 


ss 
<fet 


is the maximum velocity at y = 0, and U, is a constant 
velocity. The basic flow stream function is specified to be 
constant at the northern and southern boundaries (y = + D). 


Equation (2.10) can be integrated to give 


w(x,y) = U(x) d(x){tanh(ay) + tanh (ary) } 


- Uj(ytD) + H(-d) , (2.11) 
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where 


b(D) + 20.0 - #(-D) 
U(x) = > oer coth (Ty? . (2.12) 


Therefore, if d{x) varies slowly in x , so does the basic 


flow. The x-variation for the characteristic length scale is 





given by 
X-X ) 
850 km + 350 km{cos[27 —2 }} , 
d = L 
1200 km 
BOX eae ee 

O-=—= = Gg 
ame Se Sr : 
X > X) : (2 ilk } 


Here, A is the longitude where the x-variation of the cosine 


~ 


function starts and kt is the wavelength of this variation. 


TS 


ae 





Til. FINITE DIFFERENCE EQUATIONS 


Equation (2.7) is solved by finite differences. Centered 
time differences are used except for the friction term which. 
iS approximated with a forward step. The Jacobians are com- 
puted with the form which was developed by Arakawa (1966) for 
non-linear conservation. 

The tendency of the disturbance stream function is de- 


fined by 


(eo. ) 


+ 
— | 
vd 
qa 
li 
oe | 
Qo] a? 
ct] 
_— 
. ct 
Ld 


where superscript t is the current time and subscripts i 
and j refer to the x and y_ grid points, respectively. 
In the following finite difference equations, At({l hour) is 
the time step. Ax(375 km) and Ay(125 km) are the x and 
y grid point intervals, respectively. Equation (2.7) is 


written using the Leapfrog finite difference scheme: 


Coe a t Tce t ees 
PT DGG - Dig oe 


\ c c 
ip ie Miele oe ec ent 
ZAX acy 
(3.2) 
(teat _ |, t-At 
Ving = YElg ot 7G gAt (233) 
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The finite difference form for the Laplacian of the tendency 


is given by 





oat Tie i iat Ries inl ed Geel 
oJ 3 Ay? 
(3.4) 
and for the Jacobian terms it iS given by: 
= _attliae, aT ai, ne = c 
- (W ) ( 7G 
Yai jt “Vi me ae Meo cee : 
@ t 1 t t 
Ni, 5 7 Toxby ier, Sb41,g41 7 Fi41, 5-1) _ (3.6) 
i al rt 
We Se ee re ee 
ay ee st 
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+f it mF 
cs es ~ cH -1,j-1)! ; 
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C - 1 ac = 2 
aT oe STS gee - Vil ge1? 37) 
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+ C; 


ol 
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i-1,5 G-1, 541 esa 
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where 


TF jit vb is expanded in the same manner as for (3.8). 
This Jacobian conserves both kinetic energy and enstrophy in 
the None Ti Wear formulations (Arakawa, 1966). 

A forward time step is used for the initial time step and 
every 30 time steps thereafter to damp the internal gravity 
waves produced by the leapfrog method. A Matsuno finite dif- 
ference scheme is tested in place of the forward time step 
and the results show that the Matsuno scheme provides further 
damping, but that the change is not significant. 

The boundary conditions employed by Tupaz et al. (1978) 
relieved the problem of reflection of incident waves from the 
boundary back into the interior region (Pearson, 1974). An 
Euler-backward finite difference scheme was emploved every 
time step to provide additional damping of the short waves. 

The frictional terms are evaluated at the previous time 
step in order to insure linear computation stability (Haltiner, 
1971). 

The Poisson Eq. (3.1) is solved for the tendency using the 
POTSS subroutine which is a cyclic reduction algorithm for 
Solving block tridiagonal systems of arbitrary dimensions with 
Dirichlet boundary conditions (Swarztrauber and Sweet, 1978). 

The subroutine POISS solves the linear system of equations 


of the form: 


18 





(iy a = Y 


Nokes feats goa 141d i .d-1 o*i.5 1 ie i,j 


Pa los) T Vs 


(329) 


The east-west boundary conditions for all quantities are 


represented by 


Nig 7 MG 


Xo 4] j = 0, Jeti (3.10) 
It should be noted that Eq. (3.9) is not in the correct form 
which satisfies the model. Therefore, the model must be 
flipped prior to entering the solver and then flipped back 


upon exiting with the solved eigenvalue-solutions. 


es, 





IV. LOCAL STABILITY OF BASIC FLOW 


A parallel flow model of Williams et al. (1971) is used 
to determine the local growth rate of the variable mean flow. 
This is taken from Tupaz et al. (1978) and then transposed 
to fit the basic flow model used for this study in order to 
gain some insight on the stability characteristics of the 
mean flow. This model is hereateer referenced as the par- 
allel flow model. If d3W/dx is set equal to zero in (2.5) 


the governing equation of this model is 


fe) = 3 2 ’ a ,) ay as ‘ 
[a5 + U ay + De] Vp + [8 oe | GO (4.1) 
u(y) = -U sech® (4) ale 


where U 1s a specified constant that scales the magnitude 
of the central velocity of the Bickley jet (y=0). The 
characteristic length d is also a specified constant. 
Equation (4.1) is solved with the initial value technique 
used by Williams et al. (1971). Assuming that all perturba- 
tions are periodic in x, Eq. (4.1) is Fourier transformed in 
x and the resulting equations are integrated until the solu- 
tion grows.with a constant exponential rate. This approach 
Gives the growth rate, phase speed and wave number structure 
of the most unstable meridional mode as a function of the 


ZOnal wave number k 
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In general, these equations have a set of normal mode 
solutions and a continuous spectrum of solutions (Case, 1960; 
Pediosky, 1964 and Yanai and Nitta, 1968). However, only 
the normal mode solutions can give significant growth and 
the most unstable mode will dominate after a sufficient 
period of time. 

Figure 1 shows the growth rate (n) corresponding to the 
most unstable wavelengths (L) as a function of x based on 
the parallel flow model. The parameters for the mean flow 
are given in Section V. The largest growth rate occurs at 
x=0 where the jet achieves its maximum central velocity 
fu = - 30 m sac]: In the inflow and outflow region for 
x < =- 11,700 km and x > 11,100 km, respectively, the growth 
rate becomes negative and approaches the linear frictional 
damping rate. The most unstable wavelengths range from 
3650 km at x=0 to 4600 km in the inflow and outflow regions. 
This variation follows the variation in the y-scale of jet 
d(x) which is given by (2.16) (Tupaz et al., 1978). 

Eigensolutions of the most unstable discrete mode for 
two jet profiles are shown in Fig. 2. Only the lower half 
of the y domain is shown because the solutions are symmetric 
about y=0. Both solutions show barotropic growth as the 
phase tilt is opposite to the shear in U. The maximum tilt 
Occurs at x=0 where the growth rate is the largest. In 
addition, the eigensolution amplitude is a maximum at yv=0 for 


x=0, but at x = 9314 km the maximum has shifted to y = + 800 
km. 


2] 
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V. NUMERICAL RESULTS FOR EXPERIMENT I 


The numerical investigation of the variable mean wind 
model comprises two experiments. Experiment I is the princi- 
pal one and it employs the same friction that was used by 
Tupaz et al. (1978). A friction coefficient of Ds = 0.15 
x 107? sec’! is used which is equivalent to an e-folding 
decay time of approximately 7.7 days. The second experiment 


is the frictionless case and will be discussed in Section VI. 


The following parameters are used in Experiment 1. 


2D = 4000 km, X, = xX = 40,125 km, U(0,0) = -30 m 7 ae 
U = 0; L = 36,000 km, Ax = 375 km, Ay = 125 km , 
At = 1 hr, X> = + 18,000 km 


Figures 3 and 4 show the stream function wW(x,y) and the 
zonal velocity u(x,y) which are derived from the above 
Parameters. The model has 108 grid intervals im the x-direc- 
tion and 32 intervals in the y-direction. Tupaz et al. (1978) 
used the same domain, but his mean flow was not the same at 

xX = xe and X 5 ie However in this study, the scale of 
variation of the jet is reduced to insure that the basic 

zonal velocity would be periodic in the inflow and outflow 
regions. Equation (2.13) specifies that the jet scale, d(x), 
varies between 500 km at u(x,0) and 1200 km at u(x,0)_. 


max tate 9 a 
The initial disturbance field is of the form: 
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A wave number of 10 is chosen as a result of the average wave- 
length which was obtained by Tupaz et al. (1978). The initial 
conditions and the jet flow are both symmetric about y=0 and 
this symmetry is maintained throughout the integration. 
Therefore, only the lower half of the y plane will be depicted 
in the results of both experiments. 

The integration of the numerical model is extended for a 
period of 170 days. A close examination of wave structure, 
phase tilt, period, wavelength, and spatial growth is made 
for the segments 80-87 days, 125-133 days and 150-157 days. 
Overall growth occurs during the entire period and is expected 
Since the waves travel through the region and then reappear 
at the input boundary. It was expected that the overall 
growth would be exponential in behavior. Figure 5 shows the 
root mean square for the entire domain of w' plotted on a 
log scale vs time. If the growth were truly exponential, the 
curve in Fig. 5 would be linear. But in fact the curve is 
quasi-linear with departures from the linear behavior occur- 
ring periodically about every 56 days. Figure 5 shows the 
first departure occurs near t=21 days and continues for 28 
days. It should also be noted that the growth rate is nega- 
tive from t=38 days to t=47 days. This departure from the 
linear appears again at t=/75 days for 25 days and again at 


t=135 days for 26 days. 
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The w' field at t=30 days is shown in Fig. 6. An entire 
train of barotropic waves actually exists throughout the length 
of the channel, but the waves upstream of x = -3500 km and in 
the outflow region are not shown.because of their relatively 
small amplitude. This is the same type of wave behavior ob- 
served by Tupaz et al. (1978). 

The phase is calculated for the segments mentioned earlier. 
This is accomplished by observing each grid point in the wp' 
field and recording the first three times it passes through 
Zero in either direction. A seven day period (P) is used in 
order to insure that a full period is recorded. Tupaz et al. 


(1978) observed a periodicity of 3.25 days. The actual phase 


for each grid point is computed using the following formula: 


0 
360 (N-N.) 


an aia ( 5a } 


where N represents the first time it passed through zero and 
NG is the starting point in time of the sign calculations 

for each time segment. Although the wave continues to grow 
during the time segment, the period should not change appre- 
Ciably. Figure 7 represents an analyzed section of the phase 
tilt for the segment 80-87 days. We note that all the waves 
Shown tilt oppositely to the mean wind shear, which is neces- 
sary for barotropic instability. The maximum tilt is observed 
downstream of x=0 where the jet velocity is maximum. 


The periods at y=0 and y=+/750 km for all three time seg- 


ments are shown in Figs. 8-10. Although the periodicity is 
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not constant, it is about 77 hours (3.2 days). This is also 
representative of the entire channel. Deviations near the 
inflow/outflow boundaries could be attributed to the small 
amplitudes in this region. Figures 10a and b show an increase 
in period at x = -6500 km to 93 hours (3.87 days) and 90 hours 
(3.75 days), respectively, and a decrease at x = 6500 km to 

39 hours (1.65 days) for a and b, respectively. Figure 8 
shows only an increase at x = 6500 km to 100 hours (4.17 days) 
and 93 hours (3.87 days) for a and b, respectively. This 
abnormal behavior could be connected to the departure of the 
growth rate from the linear behavior illustrated in Fig. 5. 

Figure 11 is an overview of the phase tilt for all three 
time segments. Although they are not the same, they are 
Similar in the region from x = 3875 to x = -12000 km. They 
tend to be different in the region where the amplitude is 
small. 

The envelope of the wave packet is evaluated at y = ~/750 
km where the wave disturbance amplitude is large. The nor- 
malized w' field is used because the waves are still grow- 
ing and results are something close to a sinusoidal variation. 
This envelope, <wW'(x)>, is obtained by recording the average 
of the magnitudes of the maximum and minimum V'/Vlave 
values that occur at each longitude over a period of 7 days. 


Here, is the root mean square w' . Figure 12 shows 


Wave 
the wave packet of all three time segments for y = + 750 kn. 
Note that the maximum amplitude in Fig. 12 is near the point 


where the parallel flow growth rate drops to Zero in Fig. 1. 
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Although they are not the same, it should be noted that they 
are downstream of the max jet core. 

Tupaz et al. (1978) derived the disturbance kinetic 
energy equation to better understand the dynamics of baro- 
tropic instability. They found a number of energy terms but 
observed that the Reynolds stress term, <(u'v' du/dy)>, was 
the only dominate source for disturbance kinetic energy. This 
Supported the large amplitudes and spatial gradients of wy' 
encountered a considerable distance downstream of the jet. 

It is expected that this energy balance will also hold in this 
study. Since this is a region of tilt opposite to the shear, 
there is a conversion of mean flow kinetic energy to distur- 
bance flow kinetic energy. Figures 1 and 2 illustrate the 
large energy conversion of this tyne for the parallel flow 
model. 

Figures 13-16 show the phase angle 98 of the wave distur- 
bance as a function of y for three different time segments. 
The phase, 84? for the most unstable wavelength of the par- 
allel flow model is also included (see Section 4). tle note 
that the waves upstream of -108/75 km tilt oppositely to the 
mean wind shear, which implies barotropic instability. In 
one case, the maximum tilt is downstream of x=0, which is 
where Tupaz et al. (1978) observed the maximum tilt. Between 
x = -10875 and x = 3750,the phase tilt compares closely to 
the parallel flow model ae At the eastern boundary (Fig. 15b), 
54” compares rather well to ae but o* and 54% Show exces- 
sive tilt for the comparatively weak instability in this re- 


gion. This excessive fluctuation may partially be explained 
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by the small amplitudes which are present in this region. 


Near the outflow boundary (Fig. 16) the tilt is reversed in 
* . 


2 
energy from the disturbance back to the mean flow. This can 


S) and 05. This indicates dynamic stability or a flow of 
also be detected in Figs. 11b and c. In almost all cases, we 
observe that the tilt of the wave disturbance qualitatively 
adjusts to the local stability of the mean flow. Since the 
parallel flow model can only solve the most unstable discrete 
mode, dynamic damping is not indicated in Fig. 15 for the 
outflow region. Note that the phase tilt is unstable on the 
average. This iS because the waves show a net growth on 
moving through the region. 

The wavelengths for the numerical model are computed 


uSing the formula: 


Poel ck 
A® 


where AQ represents the change in phase over 2Ax 
Equation (5.3) can be simplified and expressed in degrees 


resulting in 


0 
eo (5.4) 


Ae° 
In addition, a smoothing process 1S applied in those areas 
where the wavelength Varies noticeably to give a more accurate 
picture of the wavelength throughout the channel. This cal- 
Culation is made by choosing a questionable point (OO) in 


the y domain of the phase field and adding 180° to it. This 
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new phase is then located by moving in both x directions. The 
distance between the new phases 1S calculated in terms of the 
number of x grid points and multiplied by Ax to give the 
correctedsor smoothed,wavelengths. Figures 17-19 show the 
variation of wavelength, L(x)’, for all three time segments 


at y = + 750 km. L(x) represents the smoothed wavelength. 


O 
Comparing the wavelengths of y=0 (not shown) and y = + 750 km, 
it is observed that the wavelengths are quite similar for 
each time period. The maximum wavelengths are within + 4500 
km of x=0 and tend to decrease toward the boundaries when 
the wavelength pattern becomes noisy. This is the same region 
that the phase tilt (Figs. 15b, 16) experiences large devia- 
tions from the parallel flow model due to the small ampli- 
tudes. Tupaz et al. (1978) observed that the wavelength at 
y = + 750 km was Jarger upstream and smaller downstream of 
X=0 than at y=0. 

Figures 20-22 contain the spatial growth, ma computed 
from the numerical model, respectively. Figure 23 is the 
smoothed version of nm. 


] 
the parallel flow model obtained by using 


3 and the spatial growth rate (m) of 


m = ene (5.4) 


where C. is the phase of the waves and n is the local 
growth rate (Tupaz et al., 1978). The value of m4 is com- 
puted directly from the wave packet envelope. In Fig. 22, 
we see that m3 has a larger maximum than m and the maximum 


is shifted slightly upstream from the jet maximum. All three 
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Curves pass through zero at the same point in which the wave 
packet envelopes are a maximum (see Fig. 12). Note again the 
noisy pattern near the eastern boundary. Tupaz et al. (1978) 
observed these same features except that the m curve was 
skewed to the left with respect to m and the jet maximum. 
This was the result of the tilt structure in the wave lagging 
Spatially behind the expected value from the local stability 
conditions. The (*) represents the maximum and minimum ob- 
served by Tupaz et al. (1978). Note the similarity of the 
maximum and minimum to a and Myo: respectively. The reason 
that m3 Shows more damping in the outflow region than the 


parallel flow model is because the latter can only give the 


eigensolutions of the most unstable mode. 
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VI. NUMERICAL RESULTS FOR EXPERIMENT II 


Experiment II is the frictionless case and employs 
the same parameters and calculations procedures that are 
used in Experiment I. However, ai ehotan the integration for 
this experiment is extended for 170 days, only the time seg- 
ments of 80-87 days and 125-133 days are included in these 
results. 

Figure 24 shows the root mean square w' for the entire 
domain. Note that the total growth is nine orders of magni- 
tude greater than that observed in Experiment I (see Fig. 5). 
In addition, the curve in Fig. 24 is more linear in appear- 
ance than Fig. 5 but still experiences the departures from 
the linear behavior during the same time interval that were 
observed when friction diag included. 

The periods at y=0 and y = + 750 km show very little 
variation and closely resembled those pictured in Figs. 7 
and 8. The predominating period which is representative for 
the entire channel is the same as Experiment I (77 hours or 
8.2 days). 

Figure 25 1s an overview of the phase tilt for the two 
time segments. The region of similarity extends from x = 
9000 km to x = -15,000 km. When the amplitudes are small, 
the behavior becomes erratic as observed in Experiment I 


(see Fig. 11). 
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The wave packet envelope, <w'(x)>, is shown in Fig. 26. 
When comparing Fig. 26 to Fig. 12, it 1s observed that the 
curves are almost identical. Since the tilt iS opposite to 
the shear, there is a conversion of mean flow kinetic energy 
to disturbance flow kinetic energy. 

Figures 27-30 show the phase angle tilt of the wave dis- 
turbance y-structure compared to the most unstable wavelength 
of the parallel flow model 80° In both cases, the maximum 
tilt was at x=0. Comparing them to Figs. 13-16 in Experiment 
I, we see a Similar structure except for Figs. 27 and 29a 
which show greater tilt between 500-750 km. Note again the 
excessive tilt of 8,* for the comparatively weak instability 


; 
Shown in Fig. 29b, and the elimination of the growth of 6,* 


2 

of y greater than -1000 km in Fig. 30 for the Figs. 15b and 
16, respectively, in Experiment I. The curve for o,™ Stresses 
the dynamic stability which can also be observed in Fig. 25b. 

The wavelength patterns for y=0 and y = + 750 are very 
Similar to those observed in Experiment [. The maximum wave- 
length occurred within + 4000 km of x=0. 

The spatial growth rates (Figs. 31-33) are near mirror 
images of Figs. 20, 21 and 23. The growth rates are a little 


larger than in Experiment I as would be expected because there 


is no friction. 


96 





"sA@p CEL-GZL (G) puke skep sg-og (e) :szuawhas auiq ayy . 
AOJ LSpOW LeItuauNnu ayy JOS<(X) ,f> SadoOLaAua yayoed aAem “TT YUualuadxy ‘92 “B14 


< (UY OOOL)X 
Sl él 6 S 3 0 or a5 GF Sele Si= oa Lom 





O°L 


O°? 


O°€ 


OV 


<(X) > 


a/ 





‘y= x 
apnzibuoy, wos ‘[==-)%9 * LPpow MOLJS LAa_LLeued ayz Jo yyHuaLtarem ajqeysun 
}SOW 9yzR puke (~--:) SAkep EEL-GZL (%°A) “(***)sfhep /8-08 (, 9) Squawbas 
owt} UOJ LBPOW [LPILuauinu ay, JO g aseYyd JO UOSLURdWOD “TY JUAaLUadxy 


+ (WH OOOL) A 





0 Oc- Ob- O9- O8- OOL- O2L- OvL- 


« (a5PUd) 46 


"42 ‘B14 


98 





0 "Wy OGLE = xX (q) pue Wy QOGzEe- = x (®) 
:Sapnzibuo, wos *(---) 6 ‘Lapow Mops pa__eued ayy Jo yyHuapanem a, qezsun 
29) “C) skep 78-08 (,'0) squauibas 

$oO uostuedwo) “T] YUuawWluadxz “gz ‘HL4 


ysow ay} pue (---*) Shep EE L-s2i ( 
awLz 4Ofs Lapow ,eoLuauNU ayy so A asey 


« (Wy QOOL) A 





0 O2- Ob- O09- O8- OOl- OZI- Obl- QO O2c- Ob- O9- O8- OOL- OZL- OFIL- 


+ (aseUd) 0 + (aseud) 40 


29 





0 “wy S/8°OL X (gq) pue wy S/8°OL- 
:sapnzibuo, uojy *(---) “9 *Lapow mol 
}SOW 394} pue (-+--) 


zLallered 3yu} 40 yabueLaaem aj.qezsun 
SAeWp CEL-GZL (x 8) *(°°*) SAPP /8-08 (xe) S}PUaUbas 
OWL} JOJ [apow [edLUsuINU dy JO g |BSeYyd Jo UOSLUeduO) 


“T] Queuwtuadx 


ees ® 
eee na 


! 





0 O¢c- Ov- O9- O8- OOL- Oct- OVL- 


+ (aseud) Ae + (aseud) 0° 


0 O2- Ob- O9- O8- OODL- O2I- 





"62 “Bly 


+ (WY OOOL) A 


OVL- 


60 






"Wy O00" Le- = 
apnyibuoy sof 6(---)%6 *japow MO{s La_peued ouzy JO yYbUua[aAemM a; qezsuN 
ysow ayy pue ¢--+) Shep EEL-GZl (,°8) Aion SAep /8-08 (, '6) Sjuowbas 
QWLY AOJ LSpOW LedLUauINU ayy JO 6 aSeYyd JO UOSLUedWOD “JJ JUaWLUadxXQ “OE “BLY 


W> OOOLe- = xX 


+ (WY QOOOL) A 





09 Ov Oe 0 Oc- Ov- O9- O8- 


+ (aseud) 08 


6] 








“Shep /8-Gg quewbas ewly wos &y'(x)u ‘azeu uymMoub [elzeds “TJ yUusawluadxy “te “Bly 


« (UW OOOL) X 


GI ral 6 9 3 QO €- 0-2 6- gt= Gi- 8i- [2- 



























62 





SL 





>sAep 
FEL-GZL Quawbas owl sos “i C(x)u Sazeu yyMoub LeizedsS “yy] Juawivadxy “ze “bly 


+ (Wy OOOL) X 
ZL 6 ) € 0 €- 9-. 6- Z2t- SlL- Ci-~ 12- 
Z°L- 


O°L- 


63 





"(S8/6L) “Le 28 zedn| Aq paduasgo uinwtutw pue (UNWEXeW Buy SZUaSaudau (x) 


*Lapow pedtusauinu lea) te sane (----) 


skep C€1-GZL AOS ¥°(X)M (q) pue (7**) Sep /g-0g JO x (2) pue Lapou 
MOLJ LALLeued ayz Wous S(---) (x)w Saqeu YyMoud ide A pence "Cc bid 


+ (Wy OOOL) X 
él- St- 8L- L2- 


aL. Ck 6 9 € 0 A ad 6- 





é-L- 


0°L- 


X 
Ww 01 ) w 


Ni © oe) O + OJ GT WN 
e s 
re 


° 
r= 


64 


[= 


( 


sus 





VII. SUMMARY AND CONCLUSIONS 


This study solved the model developed by Tupaz et al. 
(1978) with different boundary conditions. They examined 
the Penation of linear waves in a zonally varying easterly 
jet which was barotropically unstable in certain longitudes. 
The waves were forced with a period w on the eastern 
boundary and a radiation condition was applied on the west- 
ern boundary to allow the waves to move smoothly out of the 
region. After a period of adjustment the solution at every 
point varied with period w , but as each wave moved through 
the domain it experienced growth or decay in response to 
the local mean flow. 

In the present study, the Tupaz model was modified by 
changing the boundary conditions to cyclic continuity on the 
eastern and western ends of the channel. In addition, the 
mean was modified slightly so that it satisfied cyclic con- 
tinuity. Experiment I corresponded to Experiment I in 
Tupaz et al. (1978) and Experiment II was frictionless, but 
otherwise the same as Experiment I. 

In Experiment I, the root mean square stream function 
grew with time because as each wave moved through the channel, 
its net growth was added when it reentered the region. This 
growth was approximately exponential with the departures from 
exponential growth having a period of about 36 days. In order 


to compare the wave behavior in this study with that of 


65 








Tupaz et al. (1978), the fields were normalized to remove 

the overgrowth in time. A procedure was developed to deter- 
mine the phase and amplitude of the waves during a particular 
time segment. The following time segments were chosen for 
analysis: 80-87 days, 125-133 days, 150-157 days. The wave 
Structure features were not the Same for the three periods, 
but the features were very similar in the large amplitude 
regions. The predominant period was about 3.2 days which was 
very close to the forced period of 3.25 days used by Tupaz 

et al. (1978). The largest phase tilt was found near the 

jet maximum and in two of the segments the phase tilt reversed 
in the stable regions as was found by Tupaz et al. (1978). 
The amplitude of the wave envelope was a maximum near the 
point where the parallel flow growth rate was zero. The 
average of the maximum spatial growth rates was very close 

to the value obtained by Tupaz et al. (1978). This verified 
their result that the spatial variation in the growth rate 
could augment the maximum growth rate. The downstream wave- 
length was somewhat erratic, but it was generally a maximum 
near the jet maximum and it was a minimum where the jet was 
weak. The results for the frictionless case were similar, 
but the total growth was much larger. 

In general, the wave behavior with cyclic boundary condi- 
tions was very similar to that obtained by Tupaz et al. 
(1978). However, the solutions in this thesis were more dif- 
ficult to analyze because the wave structure in the three 
time segments was not the same. This was undoubtedly related 


to the non-exponential overall growths. Apparently as the 
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waves came out of the barotropic damping region and entered 
the jet region they did not have the structure to begin growth 
smoothly. In this region the wave phase structure was par- 
ticularly erratic. When the present model was used to deter- 
mine wave structure and behavior with other wind profiles, 
the structure should have been computed frequently and then 
the typical structure could have been obtained by averaging. 
This study should be extended to the nonlinear case, by 
including a forcing term to maintain the mean flow when there 
are no waves present. When waves are tntroduced the modifi- 


cation of the mean flow could be studied. 
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